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By returning to the underlying discrete time formalism, we relate spurious results in coherent 
state path integral calculations to the high frequency structure of their propagators. We show how 
to modify the standard expressions for thermodynamic quantities to yield correct results. These 
expressions are relevant to a broad range of physical problems, from the thermodynamics of Bose 
lattice gases to the dynamics of spin systems. 
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Path integrals convert the difficult problem of diago- 
nalizing a Hamiltonian into the potentially simpler one 
of summing over a set of all possible paths, weighted by 
the classical action [l|, [^j ■ They are particularly power- 
ful for making semiclassical approximations, where only 
a few classical paths dominate. Often the natural vari- 
ables for describing the path are conjugate. For exam- 
ple, one would like to describe a spin system in terms 
of paths on the Bloch sphere, even though the different 
components of spin do not commute [1] . Coherent states 
are often used in such cases, and can yield useful results 
[1, Here, we analyze the structure of such path inte- 
grals, demonstrating a practical scheme for eliminating 
anomalies which were first confronted in the 1980s [^-Q. 

The issues we address were most clearly described by 
Wilson and Galitsky [T^], who used two simple exam- 
ples to illustrate the anomalies. For example, they con- 
sider a path integral calculation of the partition func- 
tion of the single site Bose Hubbard model, Hss = 
■jn (n — 1) — iih, where n = a^a represents the number 
of Bosons, U parameterizes their interaction and /x is the 
chemical potential. This is a sufficiently simple prob- 
lem that one can calculate the exact partition function 
Zss, and find Zgs 7^ Zi.^. In particular, at zero temper- 
ature, the mean occupation number calculated from Z'^^ 
is (n' 



> = [[ij]] I which the exact result derived from Zss 
~ [[Ij + 5]] ■ Here [[x]] is the integer closest to x. 
We derive an algorithm for correcting the path integral 
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includes the classical path term Fq and a sum of residues 
over the poles of log[detG~^] which generally corre- 
sponds to the quadratic corrections to the free energy 
obtained from the continuous-time path integral, and 

[Gc^ljj = (i^Li^-S^ is tlic matrix of perturbation field 
propagators in frequency space. We precisely define all 
these terms below as we derive Eq. ([T]) and discuss tech- 
niques for calculating the correction terms. 

As emphasized by Wilson and Galitsky, our correc- 
tions are not related to ambiguities of operator ordering 



or geometric phases. Rather, they arise from the over- 
completeness of coherent states. 

The formulation of partition function as a path integral 
in imaginary time involves the expansion 



Z = Tre-^^ = ^(*o|e-^^|*o) 
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Here /? = 1/T is the inverse temperature. is any 

complete basis of the states, characterized by a set of 
parameters ^t, e.g. = ("i v) so that a\(n,ip)'^ = 
■^/ne*"^ in the coherent state basis of the Bose- 

Hubbard model. The sum J^^^ l*t) ^ ^ is the 
identity operator, of which we insert Nt — 1 = P/At — 
1 copies into the operator. We are now summing over 
all iVt-point paths in 'i'-space, with = ^n^. In the 
limit of small At one can approximate e~^^* « 1 — HAt 
and thus write the partition function in the form of a 
discrete time path integral Z — J ^It ^* ^ where 

the Lagrangian is 



Lt = - log [(*t I * 
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When the basis {|*)} is orthogonal, the first term in 
this expansion can be taken to be arbitrarily small, and 
one can approximate |*t+i) ~ (1 -I- Atdt) |^t), and by 
taking At —s- convert the problem into the traditional 
continuous form Q. This approximation breaks down 
when expanding in an ovcrcompletc basis, if the overlap 
between consecutive time steps remains finite for states 
that differ to a non-infinitesimal degree. 

As was previously noted 0, even in the face of this 
problem, the discrete time formulation in Eq. ([3|) remains 
valid. Our task is to develop a techniques for calculations 
using the discrete time path integrals, and to relate them 
to the more familiar continuous case. In particular we 
wish to find a correction of the form Eq. ([T]). 

To do so we follow standard procedure and charac- 
terize the states in terms of a saddle point solution 'S' 

satisfying 



S<Sft 



= 0, and a fluctuation ip^., writing 
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= ^ + if>f We then expand to quadratic order in the 

fluctuations Lt = Lo+V't-i2-«/'t+''/'t-^2A-«/'t+i+0 (|i/'t|)^ 
where the classical energy Lq and matrices 2^2,^2 A are 
independent of time. In terms of the Fourier components 



ipt, the partition function reads 



Z = Vip exp 
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where summation is over the frequencies cj„ 
n = ^^^^ . . . ^'"^ , yielding the free energy 

"detG-i 



F = Fa 



27r 



for 



(5) 



While Eq. has the same form as one finds in the 
continuous time formalism, the propagators G^j are some- 
what different, and the sum is over only a finite set of 
frequencies. Despite these differences, one can follow the 
standard procedure to replace the sum with a contour 
integral, using the identity 
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Here the last sum is over the poles w/ of / (w) inside 
the contour 7, and 7 is the complex circle defined by 
\uj\ = = ^. The notation Res[/(w),w/] refers 

to the residue of / (w) at w = w/ and here / (w) = 

slog 
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^ , these become 



A largely unimportant subtlety is that 

the singularities of / (w) arc logarithmic branch points 
rather than poles. In practice, when one calculates ther 
modynamic derivatives such as (n) : 
poles. 

In a continuous time calculation, the integral on the 
left-hand side of Eq. © vanishes as one takes At — )• 0. 
This integral involves ffuctuations of frequency Wmax = 
corresponding to the time scale separating consecu- 
tive time steps. When the basis \^t) is orthogonal these 
fluctuations are vanishingly small, but for an overcom- 
plete basis they are finite, and the contour integral docs 
not vanish. Straightforward algebra then reduces Eq. ^ 
to the expression in Eq. ([T]). 

Thus, calculation of the free energy requires only the 
form of log [dot Gj^] at its singularities oj = a;/ (At) 
and its value on the contour 7. We define ujf = 
limAt->o (At). Generically this limit is well defined, 
and by taking A^f — > 00, the set {ujf} and the behavior of 



/ (w) near w / become arbitrarily close to what is found 
using continuous time. In particular, F'-^^^ in Eq. ([1]) 
corresponds to the continuous-time path integral value 
of F. 

A clear example of this calculation is provided by the 
single-site Bose-Hubbard Hamiltonian. Using the coher- 
ent state basis and the field ijjf = {5nt,4>t), the compo- 
nents of the quadratic Lagrangian are 
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L2A = - [1 - AiAt] 
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and so 



log [det GZ^] = log [2 (1 - cos (wAt)) (1 - /iAt)] (8) 

has its singularity at w = 0, where the continuous- 
time structure result holds: this expression should 
be contrasted with the continuous time result 
log [det^^^Gji] = log[/32a;2]. As expected, the 
difference between the free energies simply comes from 
the contour integral, and F — F^^^ = — ^ and an 
irrelevant constant. 

The power of this approach is more readily apparent 
in the multisite Bose Hubbard model [ll| . Consider a D- 



dimensional cubic lattice of Ng sites with lattice constant 
flQ. There momentum is a good quantum number and 
one can consider Guj,k- The large w structure takes on 



the simple form 



log 
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log[2(l-cos(a;At))(l-f EfcAt)] (9) 



where ~ sin^ {kjao/2) — /i. By summing over 

the allowed momenta one finds simply. 
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plus a constant. This is the same n dependence as the 
single-site problem. 



Acknowledgements 

This paper is based upon work supported by the 
National Science Foundation under Grant No. PHY- 
1068165. 



[1] H. Kleinert, Path integrals in quantum mechanics, statis- 
tics, polymer physics, and financial markets (World Sci- 



entific Publishing Company Incorporated, 2009). 



3 



[2] 



[3] 



[4] 



[5] 



[6] 



[7] 



J. R. Klauder, arXiv preprint arXiv:quant- 

ph/0303034 (2003), 0303034vl, URL [8] 

|http : //arxiv . org/abs/quant-ph/0303034/ [ 

J. R. Klauder and B.-S. Skagerstam, Coherent states: 
Applications in physics and mathematical physics (World 
Scientific, Singapore, 1985). [9] 
A. Altland and B. Simons, Condensed Matter Field 
Theory ((Cambridge University Press, Cambridge, UK, 
2010), 2nd ed. 

W.-M. Zhang, R. Gilmore, and D. H. Feng, Re- [10] 
views of Modern Physics 62, 867 (1990), URL 
|http : //r mp . aps . org/abstract/RMP/v62/i4/p867_l 
K. Funahashi, T. Kashiwa, S. Nima, and S. Sakoda, [11] 



|http : //iopscience . iop . org/0022- 3719/19/11/014] 

K. Funahashi, T. Kashiwa, S. Sakoda, and 
K. Fujii, Journal of Mathematical Physics 

36, 3232 (1995), ISSN 00222488, URL 

|http : //link . aip . org/link/ JMAPAq/v36/i7/p3232/sl&Agg=doi [ 
V. Belinicher, C. Providencia, and J. da Prov- 
idencia. Journal of Physics A: Mathcmat- 
ical and General 30, 5633 (1999). URL 
|http : // iopscience . iop . org/0305-4470/30/ 16/008 
J. H. Wilson and V. Galitski, Physical Review Let- 
ters 106, 110401 (2011), ISSN 0031-9007, URL 
[htt p : //link . aps . org/doi/10 . 1 103/PhysRevLett . 106. 110401] 
Y. Yanay and E. Mueller, arXiv preprint 
arXiv: 1209.2446 14850 (2012), 1209.2446, URL 



Nuclear Physics B 453, 508 (1995), ISSN 05503213, URL 

|http : //linkinghub . elsevier . com/retrieve/pii/05503213950C|aS6j8t //arxiv . org/abs/1209 . 2446] 

M. Enz and R. Schilling, Journal of Physics 
C: Solid State Physics 19, 1765 (1986), URL 



